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C^ ■ Abstract 

In this paper, we prove that the inverse problems for 2D elasticity and for the thin 
r^^ , plate with boundary data (finite or full measurements) are equivalent. Having proved 

this equivalence, we can solve inverse problems for the plate equation with boundary 
data by solving the corresponding inverse problems for 2D elasticity, and vice versa. 
^ ' For example, wc can derive bounds on the volume fraction of the two-phase thin plate 

^j^ , from the knowledge of one pair of boundary measurements using the known result for 

; ' 2D elasticity [6 . Similarly, we give another approach to the size estimate problem for 

i-£^ , the thin plate studied by Morassi, Rosset, and Vessella 21 H]. 

1 Introduction 

^^ ■ In this note we would like to connect the inverse boundary value problem for thin plate to 

^fv , that for 2D elasticity. To begin, we first discuss the inverse boundary value problem for 

cn \ the thin plate. Let C M^ be a simply-connected bounded domain with smooth boundary 

do.. Supposed that the middle surface of the thin plate with uniform thickness h occupies 
$7. In the Kirchhoff-Love theory of thin elastic plates, the transversal displacement u 
f— s \ satisfies 

(N ■ div div(CV\) = in [7, (1.1) 






% 



where V u is the Hessian matrix of u, i.e.. 



V^-u 



«,12 M,22 



and C = {Cijki{x)),i,j, k,l = 1, 2, is a 4th order tensor satisfying that C € L°°(r2), 

(symmetry property) Cijki{x) = Cjiki{yi) = Ckujix) V x G 17 o.e., (1.2) 

and there exists 7 > such that 

(strong convexity) CA • A > 7|Ap V x G 17 a.e., (1.3) 
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for every 2x2 symmetric matrix A. Hereafter, for any fmiction n, uj denotes the derivative 
of u with respect to Xj . The Dirichlet data associated with (jl.ip is described by the pair 
{u,u^n} and the Neumann data by the pair 



{CV'^u)n ■ n = -Mn, div(CV\) • n + ((CV^u)n • t),t = (M, 



t),t, 



where n is the boundary normal, t is the unit tangent vector field along 9$7 in the positive 
orientation, «_„ = Vu ■ n, and u^t = Vm • t. The quantities M„ and Mt are known as the 
twisting moment and the bending moment applied on the boundary dQ. 
The 2D elasticity equation is given by 

div cr = in Q, (1.4) 

where a is the stress, which is related to the strain e = (Vv + (Vv) )/2 (T for transpose) 
by Hooke's law 

e = S(T. 

Here v is the displacement field and § is known as the compliance tensor. We assume that 
S € L°°{n), and that S satisfies the symmetry condition (|1.2p and the strong convexity 
condition ()1.3p (with a possibly different constant). The Dirichlet and Neumann conditions 
for ()1.4p are described by v and cm, respectively. 

It is widely known that (II. ip and (II. 4p are equivalent (see for example [5]) as long as 
Q is simply connected. As for the equivalence of boundary data, Ikehata [4] showed that 
when the full measurements (the Dirichlet-to-Neumann map or the Neumann-to-Dirichlet 
map) are allowed, knowing one set of boundary data uniquely determines the other one. 
Consequently, the two inverse boundary value problems with the full measurements are 
equivalent. Note that Ikehata's result was a uniqueness proof. We would also like to 
mention that explicit formulas for constructing the Dirichlet data {u, u^n) of the thin plate 
from the traction data cin of 2D elasticity were given in [H p. 158]. However, to the best of 
our knowledge, it does not seem that the relation between the individual boundary data 
{v} and {Mn,Mt} is known. 

It is the purpose of this paper to clarify the relation between the boundary data {v, crn} 
and {u,u^n,Mn,Mt}. In fact, we show that v on di} explicitly determines (Mn,Mt), and, 
following [3|, show that an on di} explicitly determines {u,u^n), and vice versa. As a 
consequence, the Dirichlet (resp. Neumann) boundary value problem for 2D elasticity 
equation is equivalent to the Neumann (resp. Dirichlet) boundary value problem for 
the thin plate equation, and thus the two inverse boundary value problems with finite 
measurements are equivalent. We emphasize that the results of this paper hold for tensors 
C and S which may depend on x (i.e., which have spatial variations), and which may be 
anisotropic. 

Having established the equivalence of the boundary data, we can solve some inverse 
problems for the thin plate via the corresponding results for 2D elasticity, and of course, 
vice versa. For example, we can derive bounds on the volume fraction of the 2-phase thin 
plate via the result for 2D elasticity obtained by Milton and Nguyen [6]. On the other 
hand, the estimate of the size of an inclusion for the thin plate studied by Morrassi, Rosset, 
and Vessella [Tj [8] is equivalent to the same problem for 2D elasticity, which was solved by 
Alessandrini, Morassi, and Rosset [1] (corrected in [2j). Moreover, we can also study the 



size estimate problem for 2D elasticity with certain anisotropic media through the similar 
result for the thin plate obtained by Morassi, Rosset, and Vessella in [9J. 

The paper is organized as follows. In Section 2, we prove the explicit relations between 
the boundary data of the thin plate and that of 2D elasticity. In Section 3, we discuss the 
application of this equivalence to the size estimate problem for the thin plate equation. 

2 Equivalence of boundary data for the plate and for 2D 
elasticity 

It is well known that when the domain is simply connected, the plate equation is equivalent 
to 2D elasticity equation, and vice versa. To describe the equivalence, let us define the 
rotational, self-adjoint, 4th order tensor M by 

MA = R^AR^ (2.1) 

for every 2x2 matrix A, where 



1 
-1 



Then the equivalence between ()l.ip and (II. 4j) is given by 



C = MSM (2.2) 

and 

a = R^(V2u)Rx (= ^(V^n)), (2.3) 

where in the context of 2D elasticity u is known as the Airy stress function. As a conse- 
quence of these two relations, we have 

CV\ = R^eRx. (2.4) 

Recall that the Neumann data for the thin plate equation are given by 

(CV2n)n • n = -M„, div(CV\) • n + {{£.V^u)u ■ t),t = {Mt),t. 

A comment on {Mt)^t is helpful. Since divdiv(CV^n) = 0, there exists a potential ip such 
that 

div(cvM = (V',2,-V',i)> 

and hence 

div(CV\) ■ n = W • t = V,t ondVt (2.5) 

since t = — R^n. Integrating (j2.5p along dVt from some xq G d^ and choosing an appro- 
priate V'(xo)) we obtain 

i) + (CV\)n ■t = Mt on d^. (2.6) 



2.1 Traction vs Dirichlet data 



Here, following [3| p. 158], we connect the traction data for 2D elasticity to the Dirichlet 
data of the thin plate. It follows from (|2.3p that 



Rj_crn = (V^n)t 



V-u,! • t 



(2.7) 



Thus by integrating R^crn along dQ, we recover Vn = [u^i,u^2\'^ (up to a constant) on 
di}. So u^n and u^t are recovered. By integrating u^t along dQ we recover u on dQ. On 
the other hand, from u and u nwe determine Vu on dQ and therefore an through (|2.7p . 



2.2 Displacement vs Moments 

Now we turn to the relation between the displacement of 2D elasticity and the moments 
of the thin plate. Since e = ^(Vv + (Vv)-^), we have 



R^eR_L 



4^1,2 - ^^^2,1 



4^1,2 - 5^2,1 
^1,1 



Thus we obtain 



divRieR_L 



^^2,12 - ^^'l,22 - 1^^2,12 
■"1,12 - 1^^1,21 - 5^^2,11 



R 



(ifl,2-^W2,l),l 
(iwi,2 - ^t'2,l),2 



|^'2,12 - 1^^1,22 
i^^l,12 - 1^^2,11 



RiV(-fl,2 - -V2^), 



which implies 



n • (divRleRx) = (R±n) • V(-vi,2 - -^2,1) 

= -t- V(-Vi,2--W2,l) 
= (-^2,1 - 2^1,2),t- 



It then follows from (12.51) that 



1 1 

■0,4 = n • div(CV^u) = (2^2,1 - -jVi^2),t-, 



and from (12.61) that 



Mi 



-(^2,1- 1-1,2) + (RleR±n)-t. 



(2. 



Observe that ^(^2,1 — "^1,2) can be expressed as 



2(1^2,1 -^'1,2) 



R 





(^2,1 - 1^1,2) 



4(1^2,1 -'"1,2) 



R_Ln t. 



We also have 



(R^eRxn) • t 



From (12.81) we thus find 



Rl 



^1,1 

if 1,2 + ^f2,l 



^2,2 



R n t 



Mt 



R 



^^1,1 Wi,2 
■V2,l ^'2,2. 
T 



Rin 



= -t • Ri(Vv)t = -(R^t) • (Vv)t = -n • (Vv)t. 
We also deduce from the definition of M„ that 

M„ = n • (R^eR^)n = t • et = t • (Vv)t. 



(2.9) 



(2.10) 



Formulae ()2.9p and (j2.10p show that M„ and Mj can be recovered from v,f On the other 
hand, we can recover v,( from Mt and M„, i.e., v^j = — M^n + M„t. By integrating v,f 
along 5r2, we then determine v on dQ. 



3 Applications to inverse problems 

3.1 Estimating the volume of an inclusion for the thin plate with isotropic 
phases 

In this section, we want to discuss the size estimate problem for the thin plate equation 
based on the similar problem for 2D elasticity. Now we assume that the thin plate is made 
of an isotropic medium, i.e, the fourth order tensor C is given by 



Cijki = --6(1 - iy){5ik6ji + Su5jk) + BvSijSki- 
The scalar B is called the bending stiffness and is defined by 



^^h^f E 



(3.1) 



I2\l-u^ 



where E is the Young's modulus, and v is the Poisson's coefficient. Both E and v can be 
written in terms of the Lame coefficients as follows: 



^^M2m + 3A) ^^^ ^ 



X 



/i + A ~"~' ' 2(/i + A)' 

In terms of these the strong convexity condition (|1.3p reads as 

/i > 7 and 2/i + 3A > 7. 
Now assume that the plate is made by two different materials in the sense that 
A = AiXi + A2X2, /i = mXi+^2X2, 



where 

{1 in phase j, 
otherwise. 

Suppose that u is the solution of the thin plate equation (II. ip having boundary data 
{u, u^ni Mn, Mf} with 

(CV\)n • n = -M„, div{CV^u) ■ n + {{CV'^u)n ■ t),i = (Mi),^. 

The inverse problem here is to estimate the volume fraction /i of phase 1 (or /2) using 

{u,U^n,Mn,Mt}. 

By (j2.2p . we find that the thin plate equation with elastic tensor ()3.ip is equivalent to 
2D elasticity with a compliance tensor S = (Sijki), where 

1 111. 

where 

2 2 

K = — (bulk modulus), /i = — (shear modulus). 

B{\ + v) Biyl — u) 

In other words, Hooke's law is given by 

1 1,1 1 ,^ , ,, 

where I2 is the identity matrix, and the equilibrium equation is 

div £7 = in Vt. 

Using the result in the previous section, we can determine the displacement v and the 
traction an on 5i7 from {u, n^„,M„, M^}. So to solve the size estimate problem for the 
thin plate, it suffices to solve the same problem for 2D elasticity. 

Therefore, when all moduli involved are known, the result obtained by Milton and 
Nguyen [6] will give us bounds on the volume fraction and also the attainability condi- 
tions for these bounds for the two-phase thin plate with homogeneous isotropic phases. 
There are two approaches used in [6j - the method of translation and the method of split- 
ting. These methods depend on suitable null-Lagragians, which are (possibly nonlinear) 
functionals of fields that can be expressed in terms of the boundary measurements. For 
2D elasticity, it is known that (c), (e), (deter), are null-Lagrangians, where (/) is the av- 
erage of the field /. Thus, from relations (12. 3j) and (12. 4p . we immediately get that (V^it), 
(CV^u), and (detV^u) are null-Lagrangians for the thin plate. It is interesting to point 
out that (detCV^n) is not a null-Lagrangian for the thin plate. This corresponds to the 
fact that (dete) is not a null-Lagrangian for 2D elasticity. 

We now consider the case where the two phases of the plate are themselves inhomo- 
geneous, with spatially varying moduli. Assume that the medium in, say, phase 1, is 
given, but the medium in phase 2 is unknown. We also want to find upper and lower 
bounds on the volume of the unknown inclusion (phase 2) from one set of measurements 
of {u, u^n, M„,Mt}. Under suitable conditions, this problem was solved in [7] (with the 
fatness assumption) and in [8j (for a general inclusion). The method used there is based 



on some quantitative estimates of the unique continuation property for the thin plate. 
Using the equivalence of inverse problems for the thin plate and for 2D elasticity, we can 
convert the size estimate problem for the plate just described to the same problem for 2D 
elasticity, which was solved in [1| (and also in [2]). 

3.2 Estimating the volume of inclusion for 2D elasticity with anisotropic 
medium 

Now we assume that 2D elastic body is made of inhomogeneous and anisotropic medium 
with compliance tensor S = {Sij^i). Due to the symmetry property for the compliance 
tensor S, we can denote 



•5*1111 

'S'lll2 

'S'2212 
'S'1212 

5*2222 



F, 5*1122 = 5*2211 = B, 

5ii2i = 5*1211 = 5*2111 = —D, 

52221 = 5*1222 = 5*2122 = —C, 
5l221 = 5*2112 = 5*2121 = E, 

A. 



Suppose that the elastic body consists of two different media, i.e., S = §0X1+^X2 • Likewise, 
the stress tensor a satisfies 

divo" = and e = Sa in 0,. 

Now assume that So is given and S is unknown. We are now interested in estimating the 
size of X2 by one boundary measurement {v, an}. 

In view of the relation (j2.2p . the corresponding elastic tensor C = C0X1 + CX2 = {Cijki) 
of the thin plate is given as 



Ciiii 
C1122 
C1112 
C2212 
C1212 
C2222 



5*2222 = A-, 
52211 = B, 
— 5*2212 = C, 

—5*1112 = D, 
5l212 = -E*, 
5iiii = F. 



So the problem above is equivalent to the size estimate problem for the thin plate with 
elastic tensor C using {n, n^„.M„, Mt}. This problem has been studied by Morassi, Rosset, 
and Vessella [^. A key ingredient in the proof of [^ is the three-ball inequality for the 
plate equation with elastic tensor Cq. The three-ball inequality is proved under a so- 
called Dichotomy condition for Cq, which permits us to decompose the fourth order elliptic 
operator associated with the thin plate equation into a product of two second order elliptic 



operators. We now write 


i this Dichotomy condition in terms of §o = 


= (5°h)- We define 




f cO _ 77 qO — qO - 

'-'1111 — -^0) '-'1122 — '-'2211 - 


= Bq, 






qO _ qO _ 
'-'1112 — '-'1121 — 


'-'1211 


- 'S'2111 - --^0, 




< 


qO — cO _ 
'-'2212 — '-'2221 — 


c-O 

'-'1222 


~ '-'2122 ~ ~^0) 






'^1212 — '^1221 — 


'S'2112 


= 'S'2121 = -^0, 






,'^2222 ~ ^0) 










ao = ^0, ai = 4Co, a2 = 


2Bo + 4Eo, 


as = 4Do, 


04 = -Fo, 


and the matrix 












oo Ol 


02 


03 


04 0" 






ao 


ai 


02 


03 04 









ao 


ai 


O2 ^3 O4 




7W(x) = 


4ao 3ai 


202 


03 









4ao 


3ai 


202 


03 









4ao 


3ai 


2a2 03 








. 





4ao 


3ai 2a2 03. 





(Here we use similar notations as in |5 
1 



Then the Dichotomy condition is defined as 
1 



ao 



|detA^(x)|>0 or — |detA^(x)| = V x e f7. 



ao 



(3.2) 



Assume that (13. 2p holds and Sq € C^'^{^) (equivalently, C*^ G C^'^{Q)). Then one can es- 
timate the size of X2 by one boundary measurement {v, an} under the fatness assumption 
(see [9l Theorem 3.2] for details). 
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